Mathematical modeling of magnetostrictive nanowires for 
sensors/devices with application to Galfenol 

Krishnan Shankar 
^^ Aerospace Engineering and Mechanics, University of Minnesota 

o 

^ Abstract 

> 

Q Magnetostrictive wires of diameter in the nanometer scale have been proposed for ap- 

^ ^ pHcation as acoustic sensors [Downey et al., 2008) , |Yang et al., 20061 . The sensing mech- 

,—1 anism is expected to operate in the bending regime. In this work we derive a variational 

^_^ theory for the bending of magnetostrictive nanowires starting from a full 3-dimensional 

* 7^ continuum theory of magnetostriction. We recover a theory which looks like a typical 

^ Euler-BernouUi bending model but includes an extra term contributed by the magnetic 

'T~| part of the energy. The solution of this variational theory for an important, newly devel- 

^ oped magnetostricitve alloy called Galfenol (cf |Clark et al., 20001) is compared with the 

. result of experiments on actual nanowires (cf jDowney, 2008J ) which shows agreement. 

g 

I 1 Introduction 

Q Magnetostrictive solids are those in which reversible elastic deformations are caused by 

I ^1 changes in the magnetization. These materials have a coupling of ferromagnetic energies 

^ with elastic energies. Typically magnetostriction is a small effect in the range of 20-200 ppm 

J> for commonly occurring ferromagnetic materials like Fe, Co and Ni alloys. In the 1970's 

^0 giant magnetostrictive alloys like Tbo,3Dyo.7Fe2 were developed. This alloy called Terfenol 

f — ■ has high magnetostriction of the order ~ 2000 ppm, but is very brittle, and has low tensile 

strength of the order ~ 100 MPa. For this reason in most sensor/actuator applications it is 
used under compressive strain. Recent research by Clark et al. [Clark et al., 2000) has led 
to the development of a new alloy called Galfenol with formula Feioo-xGcix where x ranges 
from 10% — 30%. These alloys have relatively high magnetostriction ~ 400 ppm and high 
tensile strengths ~ 400MPa. 

In recent years a lot of new experimental techniques have been developed to manufacture 
ferromagnetic wires of nanometer diameter such as electron-beam lithography, step growth 
electro-deposition, and template-assisted electro-deposition. A possible application of these 
nanosize wires is in making acoustic sensors. The inspiration for this application comes from 
the structure of the human ear. The inner ear has fine cilia like hair whose response to 
impinging acoustic waves is transmitted by the nervous system to the brain. Such biologi- 
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cally inspired devices have been proposed to detect acoustic, fluid flow and tactile inputs (cf. 
Pfang et al., 2006) ). One possible arrangement of galfenol nanowires is in the form of an array 
depicted in Fig[l]. Here impinging acoustic waves are expected to change the magnetization 
of the wire array by inducing bending deformation. 
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Figure 1: Proposed model device using nanowires of Galfenol 



The models of a vibrating string and the bending of a beam are important models in 
elasticity which are known to approximate the full 3-D behavior of a deformable body in 
the linear strain regime. Starting in the 80's rigorous mathematical methods based on the 
theory of F-convergence were used to justify these 1-D models as the correct approxima- 
tion of 3-D elasticity, loosely speaking under asymptotic conditions as the diameter of the 
3-D body approaches zero. The basic references for these results are jAcerbi et al., 1991) and 
[Anzellotti et al., 1994 1, while reference for T -convergence can be found in jBraides, 2002| . 

Meanwhile in the micromagnetics literature there has been extensive use of F-convergence 
based methods to derive reduced dimension models for ferromagnetic thin films. The earliest 
results in this direction are [ |Gioia and James, 1997| and jCarbou, 2001[ . Since our nanowires 
are expected to be used for the proposed sensor application in the bending deformation regime, 
the main goal of this paper is to combine the ideas of the references cited above from the elas- 
ticity and micromagnetics literature to derive similar asymptotic models for magnetostric- 
tive nanowires in bending. The nanowires we are modeling have diameters in the 10-lOOnm 
range with lengths in the range 2-5/im. We will show that the bending behavior of a mag- 
netostrictive nanowire resembles the classical Euler-Bernoulli bending model with an extra 
term which comes from the magnetic part of the energy. 

§|2]gives a brief review of the continuum theory of magnetostriction and defines the classi- 



cal energy S{m, w ) as a function of the magnetization-deformation pair (m , u ). The section § 
|3]gives a simple heuristic argument to show the various scales of elastic and magnetic energy 
relevant to the final result. In §|4]we start with the energy &(m,u) defined on a wire of diam- 
eter £ and on rescaling the wire to have unit diameter, recover a new energy ^^{m, u) which 
equals the energy §{m,u) per unit wire cross-sectional area, and depends on a rescaled 
magnetization-deformation pair (/re, u) now defined on the wire with unit diameter. Start- 
ing with minimizers {m^ ,u'^) of the energy J^^(m,u) in §[5] we derive the first variational 
limit problem which physically represents the magnetoelastic equivalent to the elastic theory 
of an extensible string. § [6] gives the next order correction to the first variational problem 
which only involves magnetic terms. § [7] gives the following order variational problem which 
is the main result of this paper and describes the bending behavior of the magnetostrictive 
nanowires. Here we show that we can extract a deformation w^ (cf ( |7.11| l] from the energy 
minimizing pair (m'^,u'^) which itself minimizes an energy J'^ (cf ( |7.19^ ] where J^g i^ ^^ 
energy which resembles the classical Euler-Bernoulli bending energy with some correction 
terms depending on the magnetization. The method of proof involves the idea of convergence 
of minimizers, and we do not use the more abstract F-convergence method. The Appendix \A\ 
treats the magnetostatic energy separately. 

Basic notation: a,l5,Y,--- are scalars; a,u,m,--- denote vectors in IR ; A,B,E,--- are 
tensors in IR and S c IR"' represents the surface of the unit ball in K . Components of 
any vector m are denoted by either mi,m2,ra^ or mx,ray,my. For any matrix A, A de- 
notes the transpose of the matrix. We use standard function space notation of L (n,IR ), 
i?^(IR^,K^), Hl(CL,W^); for details refer | |Adams and Fournier, 2009) . By Young's inequality 
we mean 2a6 < 5~ a + 56 for IR35>0, a variation of the classical Young's inequality. 

2 Micromagnetics 

The initial model for ferromagnetic solids was proposed in [Landau and Lifshitz, 19351 
where they also derived a model for magnetization dynamics. The continuum theory of fer- 
romagnetic materials was developed in the work of Brown [Brown, 1963) which was sub- 
sequently expanded to a theory for magnetostriction in [Brown, 1966) , where a variational 
model for magnetostriction with small strain is developed. We give a brief presentation of 
Brown's work relevant to magnetostriction in this section. 

Let Qf be a smooth bounded reference configuration in IR depending on a parameter e. In 
the following sections we fill specify this dependence. Let m(y) be the magnetization vector 
at a point y £ H^. Below the Curie temperature, the magnetization is constrained to have 



constant euclidean norm i.e., 

\m(y)\ = ms a.e. yeOj.. 

For a bounded domain, this constraint implies m e L^CQ^, rngS ), V 1 < p < oo. We extend 
m by outside Q^ whenever necessary and denote it by mxnc - ni(y)xQ,^(y) which as a re- 
sult gives rnxo., e L^iW^, K^), V 1 < p < oo. We denote by S e //^(n^, K^) the displacement 
map. The infinitesimal strain corresponding to u (y) is, ( V^ is gradient wrt. y) 

E[u](y)='^[\iyu(y) + Siyu(yf]. (2.1) 

Interaction of the magnetization with the crystalline structure of a magnetic solid gener- 
ates an interaction energy modeled by a function, (p : rngS — ' [0, oo). This energy has a finite 
number of wells (say A^) along a set of constant magnetization vectors {m } e m-sS^ where 
the index k £ {l,2,---A^} and on which without loss of generality we can set (p[m ) = 0. 
The anisotropy energy thus becomes. 



Eanis= f (p[m(y))dy. 



For cubic materials ip{m) = —j(rn\rn\ + rn\in\ + in\rn\) ^ — %(m\m\m\), which along with 
the constraint \m\ = ms gives that < (p{m) < Ki. Thus 



= / (p{m{y)) 
Jn 



< Eanis = ^{miy}) dy<Ki\n,\. (2.2) 



The exchange energy penalizes variations in the magnetization in a body and thus tends to 
prefer constant magnetizations. It is modeled as follows. 



Ee.c=d[ \yyinfdy. 



Here d is called the exchange constant. Magnetized bodies generate a magnetic self field in 
all of IR . This field hfj^iy) is given by the following equation, 

yy.[-yy^'(y) + Anm(y)] = Q VyelR^ 

h'f^(y) = -yy^'(y), 

[IV^^'^-rel] = [\-h'fi,-nW=Anm-n ondQ.^. 

[| • |] represents the jump of a quantity across any oriented surface with unit normal re. The 
demagnetization energy is generated by the interaction of the magnetization m with hfj^ and 
equals 

Eaemag(m) = ^l\hUy)tdy = -\l h'ff,(y)-m(y)dy. (2.3) 

on Jr3 z jfi 



A standard upper and lower bound for Edemag is given by 

0< Edemag(m)=—- j \h 



{y)\^dy <l[ \m{y)fdy = l\n,\ml 



(2.4) 



since |/n| = rag- The energy of interaction between an external applied field ha £L (Q,IR ) 
and the magnetization over the body is modeled by the following, 

Eappim) = - \ ha(y)-m(y)dy. 
which along with Holder's inequality gives 

-K2<Eapp{m)<K2, -^^2 = ||^a||i2(n^)||»n||^2(a^). (2.5) 

The elastic energy for the magnetoelastic solid for small strains is given by, 

Eei= [ -{E[u]-Es{m)):C[E[u]-E,im)]dy. 

In this paper by an abuse of notation, we write the above integrand as 

{E[u] -E,(m)) : C[E[u]-Es(m)] =C[E[u] -E,(m)f. 

C is a positive definite fourth order tensor. Here Es(rn) is the spontaneous strain due to 
magnetization 



m ^Es{m)eM. 



3x3 



where M^^^ denotes the set of symmetric matrices of 3 x 3 dimension. For cubic materials 



"■ sym 
it's form is 



Esim) ■■ 



2m? 



^loC^i Aiiimim2 Aiiimima 
Aiiimi7n2 A100TO2 Aiiim2'^ 



i-ioo 



Aiiimims Aiiim2m3 A100TO3 
where / is the Identity matrix in IR . The form for Es(m) and |in | = nis gives 
\Es{m)\ <K3, where \Es{m)\^ = tr[Es{mfE,(m)). 



(2.6) 



|2 . 



The norm defined used above is the standard Frobenius norm, i.e. for any matrix M, \M\ is 
defined as the trace of M M. Also C being symmetric positive definite 4 order tensor gives 
for some 7, F > 0, 

r2 I -.^ ln\^i/T^'^ ■— ,-l!i,f'2| 



rlM^I >C[M] >y\M^ 



^MeM' 



3x3 

sym' 



(2.7) 



In addition to these, energy due to external force acting on the body in the form of body force or 
surface traction is included in the general energy. However since these terms are lower order 
in deformation u, they do not affect the final form of the limit problem. For our investigation 
in this paper, we neglect this term to reduce the length of the computation. Thus the full 
energy functional for magnetostriction is, 

6 {jn,u) = Eexc +Eanis +Eapp + E el + E demag 

= [ \d\vymf + (pim)-ha-m + ^[E[u]-E,(m)f\dy+^( \h'infdy. 
JCe ^ z I on jr3 

(2.8) 

3 Heuristic Scaling of energy 



In § |3.1| and § |3.2| we start with a cylindrical domain with radius £ and length 1. We 
then show how both an isotropic linear elastic energy and the magnetostatic energy defined 
in equation ( |2.3^ scale with respect to £. The scaling of the linear elastic energy has been 
know for long in the engineering literature, but a rigorous derivation starting from a three- 
dimensional linear elastic theory is relatively recent. 
3.1 Linear Elastic Energy 

Let = {(yi,y2) £ S^CO) , y^ £ (0,1)} be a cylindrical domain of radius £ centered at the 
origin and length 1 with axis aligned along the ys axis. Let Y be the Young's Modulus, 
A = 7T£ is the cross-sectional area, and I = ?£ be the second moment of area of the cross 
section. Let (ui, 1*2,^3) be the displacements in (yi,y2,y3) directions. From the engineering 
literature we know that the extensional energy of a rod along its axis is given as 

( YA\d3U3\^dy3 = Y7l£^f \d3U3\^dy3-^0{£\ (3.1) 

Jo Jo 

where U3 is the extension of the rod along its axis. From the Euler-Bernoulli model for a beam 
bending in the direction of the Xi axis, the bending energy is 

[ Yl\d33Uifdy3 = Y'^[ \d33Uifdy3-^0{£^). (3.2) 

Jo 4 Jo 

The different scaling of the two energies with respect to £ suggests to us that a linear elastic 
isotropic energy of the form 

W'iu) = j { fi\E{uf + - 1 tv{E{u)) f }dy (3.3) 



should factor into terms which are of different orders in powers of £. Using F-convergence 
this factorization into orders of powers off has been proven in [ |Anzellotti et al., 1994| . They 
have shown that, 

W^iu) = e^Wiiu^) + £^W2{ui,U2, U4) + higher order terms (3.4) 

where uiya) = (ui, U2, uaXya) and uiy^) = -i— -— | u(y) dyidy2 is the averaged cross- 

\Be(0)\ Jb,(0) 

sectional displacement and U4(y3)=-s / {y2ii'i—yiU2]dyidy2 gives the torsional 

£2 15^(0)1 JB.iOr 

component. 

3.2 Mag^etostatic energy 

For an ellipsoidal body it is well known cf | |Maxwell, 1873 1 that the demagnetization field hff^ 
and the corresponding demagnetization Edemag for a constant magnetization m are, 

hjf^ = -AnDm, Edemag = 2?! [ Volume of body) x Dm ■ in (3.5) 

where D £ IR is called demagnetization tensor. D is independent of position y, and has 
trace 1. For non-ellipsoidal bodies supporting a constant magnetization rn, it still is true that 
hfii = —AnDm. However the demagnetization tensor Z) (with trace still 1) now depends on 
position y. The magnetostatic energy is now given by Edemag = 27r( Volume of body) xDm- 
m, where D is the volumetric average of Z). For our cylindrical domain © = ^^(O) x (0, 1), D 
is a diagonal matrix with entries 

- 85- £^ A ^ ^ 1 4£ £^ . 

D33 = ^-^ + 0{e% Dn=D22 = --^ + -r+ Oie^ 

3ji 2 2 3n 4 

See [ [Joseph, 1966[ for a simple derivation of this result. The demagnetization energy for a 
constant magnetization m = (m i,m 2,ni z) is given by 



^demag — ^^ 



r2 , 1^2 



,m^ + m2 3 "^ I f~2 . ~2^ o~2K ^ r™i + ™2 ~2 



^,,,^_^3*((^5,^l)_,^.) .,: 



2 3n\^ ^ ^' -^^ 2 V 2 

= 71^ [ml + ml) £^ + £^Qi(m) + £^Q2im), (3.6) 

where Qi(m):= ^i [m\+m\)-2m'^^ and Q2('n) := ^[ [m\ + ml]-2mf\. Thus for a 
cylindrical domain with constant magnetization we can already see the presence of various 
orders of scales in the magnetostatic and elastic energy. 



4 Rescaling 

In this section we rescale the domain Q^ depending on a parameter e to a fixed domain Cl. 
The space variable in the original domain n^ is either denoted by y or z and in the rescaled 
domain by x. The gradient operators w.r.t. y and z are denoted by V^' and V^ respectively 
and gradient w.r.t. x is denoted as just V. All variables in the original domain Q^ come with 
the tilde notation, for e.g. m while variables in the rescaled domain are plain e.g. m. For 
any vector v G IR , we will write v = {vi,V2,V3) = {Vp,V3) where p = 1,2 and Vp denotes the 
planar component of i;. Analogously the gradient operator may be denoted by V = (Vp, ^3). 

Let Qf := [ y^ e w^, 3/3 e (0, 1) ] be a domain with cross-section a*^ c U^ where cDe is any 
Lipschitz domain in 2-dimensions. While the results of all the subsequent sections in this 
paper hold for any arbitrary cross-section w^, however for the sake of simplicity we set 

a)^=B^iO)<^U^ (4.1) 

a ball of radius £ in 2-D. We rescale the domain Q^ to £1 by the following one-to-one map 

xi = — X2 = — X3 = 3/3. (4.2) 

e e 

By the rescaling Q = [xp e w, X3 e (0, 1)] where U) is now a ball with unit radius in 2 dimen- 
sions. We rescale the fields fh{y), u{y^, ha(y), and hfjj^y) using the one-to-one maps 

m(x) = m(y), u(x) = u(y), ha(x) = ha(y), h''^(x) = hfn(y). (4.3) 

The map m{x) = m{y) being one-to-one means that we can invert the rescaled magneti- 
zation mix) back to the unsealed magnetization m{y). Also while the pair (^^, fh) satisfies 
Maxwell's equation on fl^, the rescaled pair [h'^,m) does not satisfy Maxwell's equation on 
Q. However unsealing the pair [h'^,Tn) to [hff^,m), solves Maxwell's equation on Q^. Hence 
the e superscript on ^^. 

The gradient operator V^' = (V^ , Sg ) operating on m(y) or u(y) correspondingly scales as, 

V^ m{y) = - Vp mix), d'^ m{y) = ^3 mix). 

Using the scaling of gradients, we rescale the strain E[u]{y) to get a new field k'^{x) as 
E[u](y) = ^[vyu(y) + yyu(y)^] 

^diuiix) ^{diU2+d2U-i]{x) l{^d-iU3+d3Ui]{x) 

^{d-iU2 + d2U-i){x) ^d2U2{x) l{^d2U3+d3U2]{x) 

,^[ld-iU3 + d3U-i]ix) ^{ld2U3 + d3U2]ix) dsusix) 

=: K^[u]{x). (4.4) 



Substituting the above transformations into equation ( |2.8| l we get 



dx 



S'^{m,u) = £^ j \ -^\Vpm\ +d\d3m\ +<p{m)-ha-m+ -[K'^[u]-Es{m)] 

Dividing above by £^ and defining a new energy J"^{m, u) := £"^<f (m, m) we get, 

^^{Tn,u)= j \-^\\/pTn\ +d\d3m\ +(p{m)-ha-m + — [K^[u]-Es{m)] \dx 

where S^{m) is defined and bounded by rescahng the standard demag bound in equation ( |2.4> 

1 r 1 r 1 

< S'{m):=— / |^^(x)|^</x < - / \m(x)fdx=-\Q.\m^,. (4.5) 

We investigate the asymptotic nature of the problem 

{&>') inf J''{m,u), j^, = \^{m,u)eH'^iQ.,msS^)xHl{Q.,U^)'^ (4.6) 

where H^iQ,,U^) = luix) e H^{n,U^) \ uixi,X2,0) = 0, V(xi,X2) e w| enforces Dirichlet 
boundary conditions at the base X3 = 0. For the subsequent sections we also use the notation 
i7j((0,l),K) defined as 

H^[(0,l),U) = \^w{x3)eH^((0,l},R^) I u;(x3 = 0) = 0}. (4.7) 

In the next section, we will start with a sequence of minimizers {m^ , M*^) of J^'^(m, u) and 
show that we can extract a subsequence whose limit relates to the minimizers of a simpler 
lower dimensional problem J"". 

5 First variational limiting problem 

Let ( '^ ) denote the cross-sectional average of any scalar/vector, i.e. for any field a{x) set 

dix3) = j- aiXp,X3)dxp. (5.1) 

For e fixed, let {m^, u^) be a minimizer of J^'^(/ra, u). We look at the behavior of^^im^, u^) 
as £ — ' 0. For that we will first show that J''^(ni^ ,u^) is bounded above and below inde- 
pendent of £. Then we will show that from the sequence {m^ ,u^), we can extract a sub- 
sequence(unrelabeled) such that (ni'^,U3^) on the subsequence converges weakly to some 



(m^jV") in an appropriate space. This convergence will be improved to strong and the limit 



im",!)") will be shown to minimize a new functional J^" in Theorem 5.1 
5.1 Boundedness of J'^im'^, u'^) 

For an upper bound on ^'^(m'^, u^) we compare its energy with a test function (m, 0) with m 
any constant vector on nisS and M = to get, 

J''{m',u')<J''{m,0)= I \(l){m) + -C[E,{m)f-ha-m]dx + S',{m) 

m^ , , 
<K4+^\n\, (5.2) 

where the anisotropy, elastic, Zeeman and magnetostatic terms are bounded using equations 
(|2[2), (|2]6|, ^J), (|2]5) and (|43}. The positivity of all terms in J'^im", u") except possibly of 
the Zeeman energy along with equation ( |2.5t gives the lower bound, 

J^%mf,u'')>-\ ha mdx>-K2. 



5.2 Weak compactness of minimizers (m^, u^) as e -* 
The upper and lower bound on J'^(m^ , u^) gives, 

if 5 > M -^|Vpm^|^ + d I d^m'f ^dx > d j Wm'fdx. (5.3) 

Then for some unrelabeled subsequence 

||Vpm''(:«:)||^2(n)<— £^ m'' ^ m" in L^(n), Vm" —Vm" inL^Q). (5.4) 

By the weak convergence oiVm'^ix) to Vm°{x) and the lower semi-continuity of norm oper- 
ator II ( • )||^2(Q) w.r.t. weak convergence we have using equation jBA) 

||Vpm''(^)||^2(n)^ liminf ||Vpin"(x)||^2(n)^ ^^™ o"^V "^^ " ^' 
which implies 

lim m'ix) = m^ix) = m°{x3) in L^(Q). (5.5) 

Strong convergence of m*^ to m" in L (H) gives convergence pointwise a.e. for a (unrela- 
belled) subsequence. The cross-sectional average of this subsequence in'^ixa) = j m^{x)dxp 

Joj 

then converges pointwise a.e. to | Tn''{x3)dXp = in°(x3). Since from Jensen's inequality we 

Jo) 
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have |/re'^(x3)| < \ni^\ = rtig, the pointwise a.e. convergence of the unrelabeled subsequence 
m'^Cxs) to Tn°(xz) gives on using L^ Dominated convergence theorem 

m^Cxa) ^ in°(x3) inL^(0,l) as e ^ 0, Vl<p<oo. (5.6) 

Also the strong convergence of m'^ to m° in L (Q) gives convergence of \m'^-\ to \m°\ in 
LHQ) because of the fact that | \m\\^ - \m1\^ f = \m\ - m°\^ \m\ + m°\^ and domination of 
\m^-f and \m"f by mf and i £ {1,2,3} denoting any of the 3 components ofm^. Then using 
the same argument as above we can derive an unrelabeled subsequence such that 

|7^(x3)-[^(x3)=|m°|2(x3) inL2(0,l)ase-0. (5.7) 

We will need equation ( |5.7^ for showing convergence of the elastic energy in Theorem |5.1[ 
Next we prove a proposition which we need for extracting weak compactness of the elastic 
terms. 

Proposition 5.1. Given m'^ e H^{D.) and i* ^ £ Hj{Q.) using eqn. ( |5.1^ , we have the following, 

II II 2 II II 2 II II 2 

H1lL2(n) = 1^1 \\^\\lHo,i) - \\"^\\LHay 
Il^3m"||i2(n) = |w| ||a3m" 11^2(0,1) ^ ll^sm" 11^2(0), 

and fori = {1,2,^} 

Proof The first result is easily seen using Jensen's inequality 

/ \m^\ dx3 = / — / m^dxp dx3< I — / Im'^l dxpdxs = — / \m'^\ dx. 

Jo Jo I \(i)\ Jtj Jo \(i)\ Jid m Jn 



To see the second result, note for i e {1,2,3} using Jensen's inequality 

2 
dX3 



I \d3mi'\^dx3= j \d3\-— j midxpl dx3= j — ^ / dsm^-dxp 

Jo Jo I MwlJo, I Jo |w|2 |J„ 

- J77\ l^sm^l^ dxp dx3 = — \d3m'ifdx. 
Jo 1^1 Jb) \(i)\ Jn 

Integrating over co and summing over i gives us the first result. Similar calculation with u"^ 

replacing m'^ gives the third result. Noting the Dirichlet Boundary conditions on w'^(x3) at 

X3 = we get using 1-D Poincare inequality over (0, 1), 

[ \ul'fdx3 <kJ \d3Ui'\^dx3 <Kei ^^ \d3u\fdx 
Jo Jo Jn |w| ' 

where Kq is the Poincare constant on (0, 1). Integrating over w gives the result. D 
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Using positive definiteness of C in ^2.7) , ( |2.6> and bounds on ^^{m'^ , u^), we get 

<- [ C[K'[u']-Esim')f + 2 [ \Esim')f < Kt. 
rJn Jn 

Combining this with the fourth result in Proposition |5.1| we have 

I \ui^fdx3<K6} \d3U^^fdx3 <K6 I \d3ulfdx<K6l \K^[u^]fdx < Kg. 
Jn Jn Jn Jn 

Thus ||w3^||^i(Q i\^oo and due to Dirichlet conditions on u^ we get Uz'^ eH}{0,l). For an 
unrelabeled subsequence we have, 

u^'(x3}^v''{x3)mLH0,l), d3U^'(x3)-d3V°(x3)inLH0,l) (5.8) 

Already from the fact that m° and v° depends only on the X3 space variable the 1-D nature 
of the limit problem becomes evident. The magnetostatic estimate in equation ( |A.20^ gives 

S'/m') -n\a)\j \ml(y3)fdy3 = Oie) + 0(£^'^) 

which implies on using strong convergence of m"^ to m° in L (0, 1) in equation ( |5.6| l, 

lim <g'^(in'^) = lim :7r|a)| | |in!,(x3)| dx3 = n\a)\ i \m°{x3)\ dxs. (5.9) 

£-0 <^ £-0 Jo ' ^ ' Jo ^ 

The magnetostatic estimate in equation ( |A.20^ and Remark |A.2| also gives 

S'^im") = 7i\a)\ I \ml{x3)fdx3 + 0{e) + 0{£^''^) = 7i\a)\[ |m°(x3)|^rfx3 (5.10) 

and thus, 

lim S'im') = lim S'Am'') = n\a)\ ( \m'i{x3)fdx3. (5.11) 

£^0 " £—0 " Jo 

5.3 Strong compactness of (m^, u'^) and variational problem 

(5.12) 



C[E]:E;EeMl;^andE33 = s 



Set fois) := min 

Note that fo defined above in ( |5.12t can be evaluated as 



fois) = cn\s\^ -2aci2\sf :=Y\sf, 



j /■o(s(X3))dX3=Y||s(X3)||^2(o_l) (5.13) 

where a = is the Poisson's ratio and Y = f cn - 2(70 12) is the Young's modulus. 



We now state the main result of this section 



12 



Theorem 5.1. There exists a subsequence (m*^, u'^) not relabeled such that m'^ — ' m° strongly 
in H^(Q.,U^), ui^ -^ v° strongly in H}{{0,1),U) and (m° ,v°) minimizes J^°(m,v} in sdo = 
kmixs), v(x3)) e HHiO, 1), m^S^) x Hj{{0, 1), U)\ where J^°{m,v) is defined as 

/■I -j^ 

^°(m,v)= I d\d3m\ +(p{Tn) + Ji\mpf + -fo[d3V-Esao{m))-ha-Tn. (5.14) 

Jo 2 

Proof. Comparing energy of J^*^ at its minimizer (m'^, w'^) with the test function (jn° ,u^) we 
get ^^(m^, u^) < J'^im", u^) which expands out as 



II 



-^\Vpm'f + d l^sm"!^ + ^(m^) -ha m' + - [K'[u']-Es{m'}f 
< [ \d\d3m°f + (p{m'')-ha-m'' + -[K'[u']-E,(m'')f]+S'Am''). 



dx + S'Am^) 



Equation ( |5.11> gives that linif^o S^(m^) = lim£_o S^(m°). Then taking hm-sup of both 
sides w.r.t. £, cancehng common terms, and noting that m^ix) — >• /re^Cxa) strongly in L^(Q), 
we can simplify the above equation to get 

d 



lim 

£ 



-0 Jn^£^ 



m'^l +d\d3m'^\ 



f 

Jn 



1 2 

dx < I d \d3Tn°\ dx. 
n 



But weak convergence ofS/m^ in equation ( [5^ implies liminf^_n / l^sm"! < I l^a/n'^ 
which combined with the lim sup condition above gives the strong convergence, 



d3m' -^ dam" in L^{n), 



-Vpin^^OinL^(Q). 



(5.15) 



Now we show strong convergence of the elastic terms. Set s^(x) and s '^(xa) as 

s^ix) := d3ul{x)-Es^^{m''), s ''(xa) = [d3ul-Es^^im'')]{x3) = dsu^ixa) - E s^^im'') 
where s^ixz) is defined using equation ( |5.1| l. Noting fois) = Y\s\^, using Jensen's inequality 

[ Md3ir3'-E:^n)dx= [ [Y\s'\^dx< [ y\[ \s'\^dxp]dx3= [ fo{s'{x))dx 
Jn Jo Jw Jo U« J Jn 

< [ C[K'[u']-Es{m')fdx (5.16) 

Jn 



where in the last step we have used the definition of fo from equation ( |5.12| l. 
By definition .EgggCm'^) = -^(imgp- |) which using |m||^ — \m'^\^ inL^(0,l) from equa- 
tion ^5.7) gives 



Es.Jmn^Es.Am") inL2(0,l). 



S33 
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The above combined with weak convergence d^u^'^ — d^v" in L (0, 1) in eqn. ( |5.8| l gives 

dsuT^'-E^^^n - dsv" -Es^^im") in LHO, 1) . 

Then noting from eqn. ( |5.13) that I fo{s{x3))dx = Y \cd\ ||s(x3)||^2(o d and weak lower semi- 
continuity of norm in L (0, 1) gives 

I fo[d3V°-Es.Jm°))dx<\iminfl foid3U^' -Es.Jm'))dx 
Jo. i^-o Jn 

< liminf f CiK^LM^^] -Es{m')fdx (5.17) 

where in the last step we use eqn. ( |5.16t . 

To get strong convergence we will show the converse inequality of equation ( |5.17t . For 
that we need to compare the energy J'^{m^,u^) with some test function based on m° and 
v°. But the lack of regularity ofv°eH (0, 1) requires a mollification procedure. Let v (xs) G 
mo, 1) and v'^ixa) - u^Us) in H\0, l)ash^ 0. Set s'^ixa) := (^au^'Cxa) -^sggCm")). Note 
iimh^os'' = (a3i;°(x3)-^s33(m°)). Define 

eE sj^-^im°)xi + eE s^2^m°)x2 - £o s'^{xz)x\ 
vlix):- £Es22iTn°)x2 + £Es^2^m°)xi-ea s^{x3)x2 ■ (5.18) 

For i;^ defined above, K^lv^fJ-Esim") is given by 

"'^ -|(axia3s'^+xi53^,,,(m'') + X2a3^si2('»'')) 
-as'' -|(o-X2 53s'' +X2 53£'s22('n°) + xi daE s-^^'^m")) 

s'' +2e(xi a3£si3(in'') + X2 53£'s23('"°)) 

where we have left out terms below the diagonal due to symmetry. A straight forward compu- 
tation gives, ( Recall Y = cn - o'Ci2 from eqn. ( |5.13t ) 

[ C[K'[vl]-Esim°)fdx= [ h{s^)dx + 0(£\ 

Then comparing energy of the test function (/n^jl^p with (m'^, u'^) gives 



K'[vl]-E,{m°)-- 



-os 



(5.19) 



Fixing h and taking lim-sup of both sides wnt. £, using strong convergence in ( |5.15^ , and 

equation ( |5.19^ the above simplifies to 

r £ 2 r C 2 

limsup I —[K'^[u^]-Es{Tn'^)\ dx < limsup I —[K^(v^j^)-Es(ni'')\ dx 






■Ms'')dx. 
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Now taking lim/i^o of L.H.S. andnoting that lim/i^os'' = idsv°{x3) - E s^^im")) gives 

limsup [ -[K'[u']-Es{m')fdx < [ - f(i[d^v° (x^) - E s^^im")) dx . (5.20) 

£-0 Ja ^ Ja^ 

Then {b.ll) and \b.20) together give along with eqn. ( |5.12| l, the strong convergence 

lim d^ui' -^ dav" in L^(0, 1), (5.21a) 

lim [ '^[K'[u']-Es{m')fdx^ [ Ifoidav" -Es,Jm''))dx. (5.21b) 

Finally its easy to see that the strong convergence from l |5.15| l and ( |5.21t together gives, 
lim J^'im', u') = \(D\J'''{m°,v''). 

Now we show that [m°,v°] minimizes J^°(in,u) in ^o- For any [m,v) with, m £ H^{{0,l),msS^) 
and V e C°°(0, 1) n {v(0) = 0}, let us define as in eqn. ( |5l8) 



V'ix) := 



eE s^.^{m)xi + eE s^2{m)x2- £(y [d^^v - E s^^{Tn))xi 

£Es22(fn)x2 + £Es^2{m)xi - £o [d^v - E s^^{m))x2 

v{xa) + 2£[Es^sim)xi + 2Es2siTn)x2) 



Then comparing energy of J''^{m'^,u^) with J''^{m,V^) we get ^'^{Tn^,u'^) < J'^{m,V^). 
Taking lim£ of the inequality we note that L.H.S converges from above lim£_o '^^im'^, u^) = 
|w|J^°(m°,i;°). Using equations ( |5.10| l and ( |5.19> with m and V^ replacing m° and i;^ in 
the respective equations it is easy to check that the R.H.S converges as lim£_o J'^{ni,V'^) = 
\(o\^°{m, v). We thus get our minimizing principle on noticing that v e C°°(0, 1) n {v(0) = 0} 
is dense in H}{0, 1) D 



5.4 Minimization of limit problem 

The minimization of J'°{ni,v) is a substantially simpler problem than the original one. One 
can see that if the applied field is a constant over the domain, the terms (p(Tn)+7i\Tnp\ —ha-m 
behaves like an "effective anisotropy". If this is minimized over constant vector m £ rngS to 
give m", then its easy to see that {m° ,E s^^{Tn'')x3) minimizes J^°{Tn, v). 

For a large class of ferromagnetic materials, the largest energy in the "effective anisotropy" 
for typical applied fields is the demagnetization term n\nip\ which finds its minimum if 
jn° is an axial magnetization (0,0, mg). In particular for our nanowires of Galfenol this is 
true. Experimentally produced nanowires of Galfenol of 30-100 nanometer diameter show 
strong alignment of magnetization along the axis in the absence of applied fields and need 

15 




Figure 2: Left: MFM scan for several Galfenol nanowires, Right: Magnified scan of single nanowire, 
scale of the wires shown in bottom of left figure (Scans courtesy of Downey [Downey, 2008| ). 

large applied fields in transverse direction to alter this state. Experimental verification of 
these results for Galfenol wires can be seen from Magnetic Force Microscopy (MFM) scans in 
Figures |2] taken from | |Downey, 20081 . 

These scans are done for wires with 100 nanometer diameter and <110> crystallographic 
orientation with no applied field. For cubic anisotropy <110> is a local minimum of the 
anisotropy energy and gives zero magnetostatic energy contribution making it a global mini- 
mum of the "effective anisotropy". The uniformity of the scan along the wire length depicts a 
uniform state of magnetization and the bright and dark spots at the two ends are interpreted 
to be the field lines due to an axial magnetization producing net positive and negative poles 
at the ends. 

With these observations in mind, for the following sections we will assume that the field 
ha is constant. This assumption simplifies the calculation in the following sections without 
effecting the main presentation of the asymptotic limiting problem. Let us then set 

Qo = f Vim°) + n\m°/-ha-m'' (5.22) 

where m" minimizes ip(m) + n\inp\'^ — ha- m in rngS^. Then {Tn°,v°) minimizes J^" where 



16 



v'':=Es^^im")x3. Set 
u''{x}:= 



eEs^^{m°)x2 + eE s^^{m°)xi 
Esss(m°)x3 + +2£[xiEs-^^(m°) + X2Es2s(m'')) 



where we have abused notation a httle as u°{x) depends on e but does not reflect that. Then 
it is easy to check that k'^[w°] = Egim") since m" is constant. Using eqn. ( |3.6t , 



[ ^[K'[u"]-Es(m'')fdx= [ lfo[d3v''-E,^^(m''))dx = 
Jn ^ Jn ^ 

\a)\ iniJ'°(m,v'') = Qo. 



(5.23a) 

(5.23b) 
(5.23c) 



6 Second order variational limit problem 

§[5] gives a rigorous derivation of the first order variational approximation ^"(ntjU) in 
the sense that for a sequence of minimizers {m^ ,u^) of ^^(m,u), lim£_o J^'^(/re'^, w*^) = 
\(j)\^°{,m° ,v°) + o{e) with {m° ,v°) minimizing J^°(in,v) in an appropriate space. Correc- 
tors to this approximation come up as higher order theories which involve an expansion of the 
0(e) term. These higher terms can be understood as an asymptotic F- series of variational 
problems in the sense of jAnzellotti and Baldo, 1993) . 

With this in mind we define J'lini'^, u'^) := e'^iJ^'^im^, u^) - Qo)- We look at the limit 
minimization problem corresponding to J'^im^ jU'^). For this we first show that ^^{m^ ,u^) 
is bounded above and below independently of £ so that its limit e — ' makes sense. We then 
show that a limit exists as £ — ' for the quantity J'^im^, u^). Note that 

J^%m',u')-Qo=\ [ -^\V„m'\^dx +\s'Am')- [ 7i\mi\^dx] + \ [ Idldam']^ 

+ (p{m')-ha-m' + n\m'p\'^ + -[K'[u']-Es{m')f}dx-Qo 



= ^' + ^' + e 



(6.1) 
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where 21*^, ^^ and C^ are the terms in the big square brackets. 
6.1 Bounds for J^lim', u') 

Since (m'^, u'^) minimizes J^'^im, u), we have J^^(m'^ , u^) < J"^{Tn°, u°) which on using equa- 
tion ( |5.23bt along with the definition ofJ'^im, u) gives us the inequahty 

^[{m', u') < ^lim", u") = - [eQiim") + e^QsCm")) < K9. (6.2) 

The lower bound for J'^im'^, u^) requires the following technical condition. See Result 8.2 in 
| |Bhattacharya and James, 199 91 and Definition 5.2 from | |Le Dret and Meunier, 2005t to see 
other contexts where such a condition is necessary to get lower bound estimates. 

Definition 6.1. We say that a minimizer {in° ,v°) of ^°(m,v) [cf.Eqn. ^b.lA) ], satisfies the 
strong second variation condition if for any (/re(x3),i'(x3)) £ ^o there exists A > such that, 

J^''{m,v)-J^''{m°,v°) = ^°im,v)-Qo>AJ { \dsm{x3)-d3m''ix3)f + \m-m''f 

+ \d3V{X3)-d3v''{X3)f}dX3. (6.3) 

provided ||/re — m°||^i.Q -^^ < KiqE and ||i'~i'°||^i(o d <Kii£ for some £ > sufficiently small 
and KiQ,Kii arbitrary constants independent of £. 

Set M'.= m'-m°. (6.4) 

Using the hypothesis that J'°(ni,v'') satisfies strong second variation condition let us show 
the following Lemma, 

Lemma 6.1. For C^ defined as in equation (JBTl), 

e > I ^ym\ul)dx-Qo > h[\\d3M'\l2^^^ + \\M' 



r£||2 



Proof For fixed Xp = (xi,X2) e w we define fJJV{x3) := m^{Xp,X3), QJ'^(x3) := u'^ixp,X3). 
Then using the strong second variation condition on J?" we get 

^''{Wl'{x3),%{x3))-J'°{m'',v°)> A {|53(93t"-in°)|V|aJt"-m°|^}dx3. 
For fixed Xp e cd, M'^{Xp,X3) = Tl^ixs) - m°{x3). Integrating above result over XpCa) gives 
( J'''im',^l)dxp-Qo= [=y°(ln^updx-Qo = A(||^3M^||i2(^)+||^1li2(n)) 

J(I) Jul 

notingQo = f J'''{m°,v°)dxp. Fromeqn. ( |5l2) , C[K''[u'']-Esimlf > foldsul-Esssim^)] 
which gives C^ > / J^^im^, u^^)dx - Qo and our final result. D 
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We use the result of Lemma |6 . 1 1 above and Proposition | A. 6| to get 
d 



> — 

c-2 



'^p"^'\\lHa)+Wd^"^'^-f^^\"^'pfdx]+A[ II ^3M^ ||i2(a) + ||M^ 11^2(0) ) 



^ ^ II Vp^1Il2(0) + f ( ll^sM^ llL2(n) + II^1Il2(0) ) + £Qi('»'') + ^'Q2(m'') -Z)i8^2 

(6.5) 
> -Ki2£ 

where Qi(m°) = |^ ( |in° 1^ - 2\mlf'j and Q2(m'') = n^[^^ - |m° 1^) as in equation j^ - 

6.2 Convergence of y^im^, w^) 
Theorem 6.1. We have the following convergence, 

lim^f(m^«^) = Ql = ^|m°|2-^|m°|2. 



.0|2 



Proof. Dividing equation ([63) by £ gives 

Taking liniinf£^o above to get, 

liminf^f(ln^w")>— [|in°|2-2|mS|2]. 
£-0 3 V ^ "^ ) 

To get the reverse inequahty we take divide eqn. ( |6.2| l by £ and then take lim sup to get, 

limsupJ?l(m^M")<limsup-(eQl + £2Q2l=Ql = ^(|m°|2-2|m^|2l. 
£-0 £-0 £ *• ' ci V '' ; 

The lim sup and lim inf inequality together gives our result. D 

Remark 6.1. In the limit we get that ^^{Tn^,u^) converges to a fixed quantity Qi(m°) de- 
pending only on m°. Qi(m°) consists of the mutual interaction of the poles generated by /re" 
on one end cj{0) of the wire domain Q. with the other end toil) giving the term -g-Trlmgl , and 
the self-interaction of the poles created by m" on the curved surface dw x (0, 1) giving the term 

-fTTlm^p. 
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7 Third variational limit problem 

As in the previous section we first define J'^im'^, u"^) '■= £'^{J''^im^, u^)-Qo-£Qi{in'')). 
We will show that J'^im^ , u^) is bounded above and below independent of £. Then we define 
w^ in ( |7.11^ and prove a weak compactness result for it. The convergence is improved to strong 
where we also define a new variational problem J'^ and show its relation with 



in Theorem 



7.1 



\d\dim'f' + q){m') 



Recalling 21*^, 25*^ and C^ from equation ( |6.1t in §[6] , we note 
J?^(m^w^)-Qo-£Ql(m°) 
= j ■^l^pTii'l'^dx] + [S'^im')-nj Im'/dx-eQi] + [ [ { 

= 2l" + ('B"-eQi(m'')) + r. 

7.1 Boundedness of J'^im'^ , u^) 

To get an upper bound on J'^im'^ , u^) we use the upper bound on J'^im'^, u^) from eqn. ( |6.2> 
and subtract Qiim") from both sides to get 

^|(ln^M^)=-[J?l^(m^M^)-£Ql] <-[^l^(m^w'')-eQl] = Qiim") < if 12. (7.1) 

To get a lower bound on J'^im^ , u^), we subtract £Qi(ni°) from the lower bound in the pre- 
vious section in equation ( |6.5^ to get 






> -ifi3£'. 



L2,n)+ 2 I ll^3M^lll2(o) + ll^'ll22(„))+^'Q2(m'')-Z)i8£^ (7.2) 



The upper bound 7.1 and lower bound in 7.2 together give with Sobolev inequality on Cl 



Ku£^ 



d 



\yr,M' 



A 



( WdaM' 



2^2 II "P"' llL2(n)^2 
- 2 il \\HHn)-^<i¥*^ WLiiQ) 



L2(n)^ll^'* 11^2(0) 

Vl<g<6 



(7.3) 



with Cq being the appropriate Sobolev constant. 
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7.2 Weak convergence of w^ 

In this subsection we will extract some energy terms from the elastic energy and define a new 



variable w^ from the extracted terms. For this we need an improvement on Lemma 6.1 . For 
that first note that using a truncated Taylor Expansion we write Egim^) as 

Esim") - Esim") = E'^im") ■ M' + -E'!.{m°)M' ■ M' + o{\M' f-) := MM') (7.4) 

where we recall m'^ = m" +M'^ and E'^{m) and E'^{m) are the 1^* and 2^^" derivatives of 
Es(m) wrt. m. Since Egim) is a polynomial function of #re, both E'^ and E'^ are bounded in 
L°° for \m\ = nis- Then using ( |7.3> we get 

A(M')<K^5\M'\+Km\M'\'^ and ||A(M^)||^2(a) s i^ive^- (7.5) 

Set u' = u'' + V. Note k'[u'] = k'Iu"] + k'UJ"]. Eqn. (|2j) and Young's inequahty gives 



4r^ 



|2 , Tl 



C[MM')]:k'[U'] <T\A(M')\\k'[U']\< \A(M')\ + -\k'[U']\ . (7.6) 

From ( [E!23l l note that k'^Lw"] -^sCm") = 0. Then equations (|2[7|, ( [532) and (|7]6} gives 






c, 

2 



1 r, ,9 4r^ , ,9 

= ^Wzv° -E.^^im')] + ^KWf A(M^) '. 



(7.7) 



Then using equation ( |7.5| l, and above result ( |7.7^ , we improve Lemma [63] to get 

e= I ^d Idam'f + (p{m')-ha-m' + 7i\m'p\'^ + -C[K'[u']-Es{m')fjdx-Q 
>j [d \d3m'\'^ + (pim')-ha-m' + 7i\m'p\'^ + -fo{dsv''-Esssim')]jdx 



,2 4r2, 



'[C/nr-^^|A(M^)|Hc/x-Qo 

r ' 






(7.8) 
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where we have used the strong second variation condition on J^" in the last step. Using this 
we revisit the lower bound equation ^7.2) using Proposition |A. 6| to give 

J?"(ln^M")-Qo-£Ql('n°) = 2t" + («B"-eQl(m°)) + £" 



The upper bound ( |7.1| l and the above equation then gives 
Set w" = (u^p - M° ,£"^(w^ - w°)) = ([/^,£"^?7|) and note 



(7.9) 



(7.10) 



kHU'} = 



hdiw' 



^(diwl + d2wl) \^2W2 \{d2wl + d^w^^) 

\(diwl + dzw\) \(d2wl + d^wl) ed^wl 



(7.11) 



Note 



Xiw'^) 



( |7.10> gives, 



|£(ic'^)| where E{iv^) is the elastic strain of field iv''. Korn's inequality in 



f |'^!(^ 'dx= [ 1^^ 'dx> [ \Eiw^fdx>a[ {\ywf + \w^\^)dx 



ifl9> 



where a(Q.) > is the Korn's constant. These results together imply for some unrelabeled 
subsequence 



E{w')-E{w°) in L\n;U^) 
Note from fOI) , 

-^ = ^Eijiw') for {i,j)e {1,2}, 



— -v°inL^(n;IR^). 



^ = -Ei3(w') for ie {1,2} 

£ £ 



(7.12a) 
(7.12b) 



which together imply after using lower semi-continuity of norm wr.t weak convergence 

\\E^Jilv'')\y^n) ^ li^o^f \\EiM)\y^a) = (7.13) 

when i e {1,2} and j e {1,2, 3}. 
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Lemma 7.1. Let the strain corresponding to a displacement field w° £ H^ (O), E{w°) be such 
that Eij(w°) = Q for ie {1,2} and j e {1, 2, 3}. Then 3y°(x^) e HUQ, 1) so that w° is given by, 

ivlix) = ivlix^) e HliO, 1) , w^iix) = w^ixz) e HUH^, D , 

wlix) = -x\ di,w\{xz) - X2 d^wlixi) + Y°(x3) (7.14) 

and H'^p, l) = {we HHO, 1) : w{x3 = 0) = dawixa = 0) = }. 

Proof Given £'ii(h;°) = diw^ix) = and ^22(1*''') = d2W2(x) = together gives wlix) = 
ai(x2,X3) and u^gW = a2tei,^3)- Ei2{iv°) = gives us, 

d2wl{x) + diW2ix) = d2ai{x2,X3) + dia2{xi,X3) = 0. 

This imphes d2ai{x2,X3) = -dia2ixi,X3) = fiixs). Thus, w^ix) = yiCxs) + X2/3(x3) and 
w^ix) = 72(^3) - xi)6(x3). Also given £13(10'') = E23{w°) = 0, we have 

d2Ei3{w°) - diE23(w°) = -{d32wl + di2wl - d3iwl - di2wl) = d3l3{X3) = 0. 

This gives us /^(xa) = -K'20 is constant. Using the Dirichlet boundary conditions at the base 
X3 = 0, we have a'^(x2,0) = 7i(0) + X2i^20 = which gives us 1^20 = 0- So w°y{x) - 7i(x3) and 
w''^x) = 72(x3). We finally have, 

Ei3{w°) = ^ diw^ix) = -d3wl{x) = -537i(x3), 
-^2,3(1*°) = ^ d2wl{x) = -d3W2{x) = -^372(^3)- 

which gives us for w^ on integrating above equations 

W^ix) = -Xi 537i(X3) - X2 3372(X3) + 7''(X3) = -Xi asif^Cxs) - X2 53"^2(^3) + 7°^), 
d3W^(x) = -Xi (333"^ite3)-^2 ^331^2(^3) + ^37°(^3)- 

Note also that w" e H^{Q,) gives 53^;° e L^{Q.). Equation ( |7.14> gives then that 533«;°(x3) e 
LHO.) for i = 1, 2 and ^373 e LHO.) and thus if^Cxs) e HH^) and 7 e H^{Q.). Note also that 
the Dirichlet boundary conditions at the base gives w^{Xp,0) - —X\ 53u;^(0) - X2 d3W°,^{!i) + 
j°{^) = which means d3w\{'d) = d3W°{0) = r°iO) = 0. D 



The displacement solution in Lemma pTl] are well known in literature as the Bernoulli-Navier 
displacements. See Theorem 4.3 in [ |Trabucho and Viano, 1996) for more details. 
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7.3 Strong convergence of J^^im^, u^) 

For the third hmit variational problem we assume that m° = {0,0, rus)- This assumption 
greatly simplifies our final limit problem while essentially describing the underlying physics. 
Refer to Remark |7.1| to see more regarding this assumption and the more general case. From 
§J^, M' = m'- m" which gives \m'f - Im^l^ = \M'f + 2m°-M' = IM^I^ + 2m°M^ = 0. 
Also note ^33(1*'^) = ^3^3 + IggCw'^) = d^v" + xls^^^^ = Eg^^im") + X33(h;'^). Then 

\M'f = -2m''-M' = -2msM'^, (7.15) 



^.33(m^) = ^(|m^3|2-^)=^.33(m'')-^ |M^|^ (7.16) 



^(K|2-^)=^.33(m'')-^'--2 
2mf ^ ^ 3 '^■^ 2mf 



and using Holder's inequality and < |7.3^ 

[ h[KUu')-E,^^{m^))= f f,[dsv''-E,^^(m^)) + fo{xls(iv')] + ^^\Mlfxls(^^) 
JQ. Jn ^"'s 

> ( fo[dsv''-E,^^{m')) + f foixU^')) -if2i£'||/33(«'')|lL2(0)- (7.17) 

Using ( |5l2) C[K'[u'']-Es{m'')] > fo[Kl^{u'')-Es^^{m'')) which gives 

e> I d \d^m'? + ip(m')-ha-m' + 2\m'/ + ]-h[Kl^(u')-E,^^{m')]-QQ 
Jn "^ A 

> [ d \d3m'\^ + cp{m')-ha-m' + 2\m'/ + lfo[d3v''-Esss(m'))-Qo 
Jn '^ z 

+ f foixl3i^'))-K2ie^xl3("'')\\LHn) 
= j ^''im\v'')dXp-Qo] +j foif33("^'))-K2i£^\\xU^')\\LHn) 
= M\\dsM'\\%in)+W'\\hin)]+ [ fo(xl3(^'))-K2i>;'\\xl3(^%mn) (^.18) 

o LI 

where we have used strong second variation condition on ^°{m'^ ,v°) in the last step. 
Define J'2{wi(x3),W2{x3),v{x3)) in function space s^2 as 

J'2(Wi,W2,r)=7: I \fo{xid33Wi(X3))+fo{x2d33W2iX3)) + fo[d3Y(X3)] dx + Q2(m'') 

AJn'- J 

(7.19) 

and ^2 ■■= {iwi{x3),W2ix3),rix3)) e ^J(0, 1) x H^iO, 1) x H^{0, 1)}. 
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Theorem 7.1. There exists a subsequence w'^ not relabeled such that w'^ — ' w° strongly in 



H-^{Q,U^). w° is given as in Lemma 7.1 and {w'^{x3),W2{x3),Y°ix3)) minimizes J'2 '-'^ -^2 and 
lim£-_o J^2^(/ra'^,M'^) = J'2iwl,W2,T'') where {w1,W2,T°) minimizes ^2^wi{x3),W2{x3),vix3)) 
in ^2- 

Proof. Because m" - (0,0, mg) we use the relevant magnetostatic estimate from equation 
( |A.23^ in the remark [AT3] following Propo sition | A. 6 1 and equation ( |7.18t to give 

J?"(ln^M")-Qo-£Ql('n°) = 2t" + («B'-eQl(m°)) + C" 

> £ { ^WpM'f + MdaM'f + A\M'f + IfoixU^')) } -K2i£'\\xk(^')h2^n) 

Jn 
- L ^^o(^33(«^')) -^2i£' II ^^^ L^fn, + ^'Q^Cm"). (7.20) 

Dividing by £^, noting e'^xls^^^^ = d^w^^ and daw^ — dsw^ in L^(Q), we get on taking 
liminf£_o, 

liminfJ^■|(m^w")> [ -fo{d3wl)+Q2im''). (7.21) 

The first term in the R.H.S comes from the fact that ( |5.12^ gives, I /'o(53U^3) = Y || 530;^ 11^2(^1 
and^aifg ^^au^g in L^(Q) which implies ||'53W^3||£2(n) -l™™f ||^3"^3llL2(n)- 
To get the limsup inequality, we compare energy of J^g ^^ its minimizer (/n^, M*^) against a 
test function {m° ,U = u° +(W p,£Wi)) where 

Wi = u;i(x3) - e^oxidzjixz) + £^-^ [xld33w'l - xldaawl + 2x-iX2d33W2] (7.22) 

W2 = w^ixs) - £^cTX2d3Yix3) + ^^ o [^2^331^2 " xldaaw^ + 2x1X2^330;^) 

W3 = W^ix) = r°(X3) - Xi53U;^ - X2d3W2 

Then K'im-Esim") = K'iu'')-Esim'') + x'(W) = x'(W) where £-'^x'(W) converges as: 
e"'lii(W^) = £"'l22(W) = -ff (53r°(x3) -xi533«^? -X2a33«;^), e"'li2(W) = 

£~^Il2(^) = 2! ~'^^1^33r°(^3)+ 2 (^1^333"^? -x|a333U;i +2xiX25333U;2)| = 0(e) 
£~^Xl3(W) = 2 { ~ 0'^2533r°(^3) + 2 (^2^333"^2 " ^1^333"^2 + 2xiX25333M;^)| = 0(e) 
e~^;t33(W) = ^3^3 = a3r°(X3)-Xia33U;i -X2a33t^2- 
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Its easy to check that C[K''{U)-Esim°)Y = fo{ed3W3) + Oi£^) gives 



J'l(m',u')<J'l(m°,U) = 



£2 



= Q2{m°)+ \ fo[diy''{xi)-xidzzwl{xi)-X2dzzivl{xz)) + 0{e). 
Jn 

Taking limsup as £ — >• we get our result. We finally need to show that {w'l,W2,T°) min- 



imizes ^"{101,102, y) in ^2- Again as in Theorem 5.1 we start of with smooth {wi,W2,y) 
satisfying the boundary conditions. We set up a displacement W exactly as in eqn. ([7]22) 
with (wi,W2,y) replacing {w'^,W2,y'')- Comparing energy ot J'2^Tn^,u'^) with the test func- 
tion {m'',u'' + (W^, W|)) we get J^^{m', u') < J?^{m'',u'' + (W^, W|)). We get our result on 
taking limit as £ — >• and noting that smooth functions (u; 1,0^2,7) satisfying the appropriate 
boundary conditions are dense in .5/2- □ 

Remark 7.1. The assumption m" = (0, 0, mg) is not necessary, but it greatly simplifies the 
form of the third variational limit problem J'^- If m" is a more general constant magneti- 
zation, then the magnetization m" will have a non-trivial corrector. In fact the third vari- 
ational limit problem then will be a more complicated problem involving both elastic and 
magnetic corrector terms. The elastic part of the problem will however still retain the Euler- 
Bernoulli type terms and the problem however will simplify to the limit problem of Theorem 



7J\ifm'' = {0,0,ms). 

We however do not present that result here, as we are more interested in nanowires made 
of Galfenol. For these wires made of Galfenol, as expressed in § \5.4\ the demagnetization term 
n\inp I is the largest term in the "effective anisotropy" (pirn) -\- n\mp | — ha ■ ni by an order of 
magnitude for typical applied fields. Thus the minimizer m" of this "effective anisotropy" is 
expected to be at or very close to (0, 0, nis). 

8 Summary and Discussion 

We have presented in this paper the derivation of simple models for nanometer diameter 
wires to be used in sensors/devices using the physics of magnetostriction. Though the starting 
point for these problems is an infinite dimensional variational problem with a non-convex 
non-linear constraint and variational energy contains terms which are non-local, using the 
method of variational convergence we have derived much simpler 1-dimensional models which 
is expected to approximate the actual physics of the starting model. The Theorems |5.1| and 



7.1 clearly set up these simpler models J'°{Tn,v) and J'S{wi,W2,v) respectively. 
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Figure 3: Bent wires of Galfenol 



The bending behavior of the nanowires is described by J'2iwi,W2,v) if m" = (0,0, mg) 
solves the first variational limit problem ^"(m, v) as is expected for Galfenol wires. The form 
of second variational limit and the third variational limit suggest that the magnetization 
remains strongly stabilized at m° and higher order theories do not add correctors to m" 
within the framework of geometrically linear theory of magnetostriction. The displacement 
solution K° corresponding to the first variational problem is however corrected due to the 
appearance of the bending energy terms in the third variational limit J'2{wi,W2,v). 

Although we have not included any external applied force in our analysis, it can be in- 
cluded with very minor changes to our presentation. The galfenol wires in bending behave 
like purely elastic beams with additional magnetic term which comes thorough the interac- 
tion of the positive and negative poles created at the two ends of the wire by the magnetization 
m" = (0, 0, rtis)- This contribution is a fixed energy at the order at which bending elastic terms 
appear. 

The strong stabilization of the magnetization is borne out by experiments where nanowires 
have been bend using an AFM tip. The Figure [3] shows the MFM scan for a galfenol wire in 
bent shape. The details of the experiment are available from | |Downey, 2008) . The MFM scan 
shows the same bright and dark spots at the two ends of the wire characteristic of axially 
magnetized wires as seen in Figure|2] The bright spot in the middle was detected to be a topo- 
logical defect. It is clear that even the large bending is unable to alter the axial magnetization, 
which can be interpreted as being equal to m". 

The bending behavior of the nanowires will be more complicated if m" ^ (0, 0, ±ins) solves 



the first variational limit problem J'''{m,v) as mentioned in Remark 7.1 This case is however 
not very important for Galfenol nanowires with the geometry that we are interested in. 

The highly nonlinear deformation of the nanowires in Figure |3] also suggests to start of 
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with a geometrically nonlinear theory for magnetostriction. For geometrically nonlinear de- 
formations however, the problem is significantly harder as the magnetic energies in the start- 
ing energy ( |2.8| l will be defined on the deformed configuration, while typically in nonlinear 
elasticity, the free energy is defined over the reference configuration. 
Recall the energy ^2(wi(x3),W2{x3),v{x3)) was defined in the previous section as 

uf2(tx)i,W2,v) = j -\^fo{xid33Wi) + fo{x2d33W2) + fo(d3r)\dx + Q2(m''). 

Note that the first and second term are exactly the bending energy that appears in classical 
Euler-Bernoulli theory. To see this note that from the definition of fo in ( |5.13t we get, 

I foiXld33Wi{X3))dx = j j Yxl\d33Wi{X3)\ = j YI j xldXp>\d33Wi{X3)\ dX3 

= I Yl22\d33W-i{X3)\ dX3 
JO 

where 1 22 is the polar moment of inertia. 

From the point of view of using Galfenol as a potential material for sensor application, the 
strong stabilization of magnetization m" = (0,0, wig) is not encouraging, as a designer would 
hope that the magnetization would change drastically from m° on imposing any bending 
deformation. Newer proposals for sensor design using Galfenol have been made which replace 
the wire array of Galfenol with an array where each wire is multi-layered with fine layers of 
magnetic Galfenol and non-magnetic Copper (cf [Park et al., 2010| ]. 
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A Magnetostatic calculations 

A.1 Introduction 

Recall in ( |4!5) we defined ^^(m) = — I |^i,(x)pdx = ^ — / \hfi,{y)^dy. In this section 
we will work in the unrescaled magnetization m and demag field ^^. We define 

m\yi):=\ m\yp,yi)dyp. (A.l) 

Note on rescaling rh'^ as in ( |4.3| l m'^ also corresponds t o the cross-sectional average defined 
in ( |5.1> i.e., m,^ - j m'^{xp,X3) dxp. Thus proposition 5.1 on m^ gives 

ll"^l£i(o,i) = II"^1Il2(o,1) + II^3'"1I12(o,i)- j^^ ll'"1lL2(n) + ll^3m^||j2(n)). (A.2) 

Let ^^e solve Maxwell equation for m'^ on Q^. We first prove the following Lemma to esti- 
mate the difference in magnetostatic energy between m*^ and m^ . 

Lemma A.I. The following inequality holds: 



1 



I. 



mHlL2(R3)-||^mHlL2(K3) -^'^I^^^'IlHO.) \^ P^^'^IlHUY 



Proof. First we recall the basic demagnetization energy bound in equation ( |2.4| >, 

-^Wm\\l2iuS)^2^\m\\l2^n^y y meL\Q.„msS\ (A.4) 

We have || in^ -ni^ IIl2(w'^) - -^i^ || ^^^^ IIl2(w'^) using Poincare inequality on a cross-section 
plane WtCya), which on integrating on 3/3 e (0, 1) gives 

Since ^gj^e satisfies Maxwell equation for m*^, by linearity (hfne —hjfie) satisfies Maxwell 
equation for ( iri^ - m^ ). Then using basic bound eqn. ( |A.4| I for Maxwell equation we have, 

■*■ lit'' r'^ l|2 ■'■||'-^£ —-Ell^ -^1 2||v7?/~£||2 

^||«m'^-«mHlL2(K3)-2ll'" " *" llL2(r!£) < ^ ^ "^"^ llL2(n,)- 

Using triangle inequality we also have. 



*mHlL2(K3)~ll^m' 



II ■^ C ■^ C II _. II ^^^If . -c II 

- ||"m'^-"mHlL2(R3)<-'>'2£||Vp/n ||jr,2(n^). 
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Jensen's inequality gives \in | < |jn | and using basic demag bound in eqn. ( |A.4| I again for 
m '^ and m'^ we have, 



Sn 



l^mHlL2(R3)+ll^^ 



Then on combining the two and rescahng m*^ to m^ we have 



l^mHlL2(R3)- 



ll^mNLa' 



"■mnlLa 



|^m'^|lL2 + ll^mnlL2 



<Z)2£||V^m^||^2(n,)- Stt ||m^||^2(n,) 



SnD2£ ^m" 



\LHa) Irp"* llL2(n)' 



by noting that ||m^||i2(n^) = f ||m^||i2(n) and IIV^m''||i2(n,) ^ IIVpm'^||i,2(n)- 



D 



Remark A.I. From ( |5.3^ we ^wow i/iai i/ie exchange energy of magnetization m^ is bounded 

as K5 > ^ II Vp/n^ 11^2(0)- -Si'ice <?J(m) = ^ll^mll^2(R3) = ^II^^IIl2(r3), we then get from 
the Lemma 

\S^{m'')-S'^{rn')\ = -^^^ ||^m'=|lL2(K3)- ll^mHlL2(R3) 



87l£^ 



<Z)2ll'n"|li2(a)||Vpm"|li2(n) 



a 



Thus the difference in magnetostatic energy between S^im'^) and S^{m^) is of order 0(e). We 

1 II "^ II 2 

will see that for the convergence arguments it is enough to estimate S'^im'^) = - — 2 ||^m''llL2(R3)- 
For miya) note V^ • m{y) = 5gm 3(3/3). It is well know that for magnetization m'^, the energy 

1 w ^^£ II 2 

- — 2 \\hfii£ ||l2(r3) can be written as a convolution of fundamental solutions with m^, 



8n 



*mMlL2 



-// 



V2'-in'^(y)V?'-in^(z) 1 
f^ly-zl ^2 

V2'-m^(y)in'^(z)-n(z) 



/ /- 

Jafif Jane 



m''(y)-n{y) m''(z)-n{z) 



e^\y-z\ 



= -( (- 



a^m 3(3/3) a^m 3(23) 



£^|y-Z| 



2 Jdn/, Jdne 



m''{y3)-niy) m''(z3)-re(z) 



e^|y-z| 



dl^mUya) m''(z3)-n(z) 



1 -.,_,, 1 



e^ly-zl 






(A.6) 



Note that ^^,^2 ^^^ '^3 respectively represent that "Bulk-Bulk", the "Boundary-Boundary" 
and the "Bulk-Boundary" terms of the magnetostatic energy. The body Q^ = a*^ x (0, 1) and 
the boundary dQ^ can be decomposed as as ^Q^ = {dcDe x (0, l)}Uw£:(3'3 = 0)U<^£(3'3 = !)• 
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A.2 Estimates of J^(to , J^ifn') , and J^iffi^) 

The magnetostatic estimates in this section are inspired by simihar estimates in other works 
hke | |Kohn and Slastikov, 2005) and [ |Carbou, 2001) . We use the following integral inequality 
in this section: for arbitrary aj^beU and q,L e U using the fact that qiq^ +L^)~^'^ < 1 we 
have 



r^ -Jl ? "-jLf ^ ? ]<A (A 7) 

i lL2 + o2l3/2 L2(L2 + ^2)l/2 „ ^2^(^2 + ^2)1/2 (L2 +a2)l/2 J " L2 ' '■' 



We also need an estimate of the following term, where we use the change of variable Wp = 
yp-Zp, dwp = dyp to get, [Recall w^ is a ball of radius e in 2-d 



dyp dzp r r dwp f f dwp 



Jai^Jaii: l^p 21^1 Jai^ Jai^-Zp lH^pl Jai^ Jo) 



«3£ I^^P 
2n rSs 



= dZp 

J(i>F Jo Jo 



Wp\ d(\Wp\)d6 



= (7re^)(27T)(3£) = 67r^e^, (A.8) 

where we have used the fact that {(De — Zp)c^ co^^ for Zp c lOe. Henceforth we drop the y 
superscript on the derivative operator 

Also note that if m'^ £ H (0, 1), Sobolev embedding gives along with (|AJ2} 



sup \m'{y3)\ < DaWm'WHHo,!) ^ —:{ ll'"'llL2(n)+ ll^3m"||^2(n))- (A.9) 



ys 



0) 



Proposition A.I. 

\Jl{m')\ < D4£{\\m'\\l2^^^ + \\d3m'\\l2^n))- 
Proof. Recalling definition of J^ from Equation ( |A.6| I and noting \yp-Zp\ < \y-z\ we have 

JneJne \y-z\ Jo Jo JoJeJlOe 



Lf:J\i£ \y Z\ Jo Jo Jw^JoJi: l^p Zp\ 

rp *pl 



Jo Jo Jcoi;JcDE\yp~Zp\ 



where we have used Holder's inequality on the term I I \d^m aiya) d^m3{z3)\dy3dz3 and 
equation ( |A.8^ in the last step. Using equation ( |A.2^ we get our result. D 
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Proposition A.2. 

\Jl(m'}\<D5£[\\m'\\l2^^^+\\d3m'\\l2^^^). 
Proof. Recalling definition of J^ from Equation ( |A.6| I, we split of J^{m'^) into 2 parts, 

J Mi; Jo JdcosJO \y — Z\ JcOeJO 

m^izs = 0)-n{z) f m^{z3 = l)-n{z) 



r r m'iz3 = 0)-n{z) C 

^JojAO) \y-z\ Ja 



ojao) \y-z\ JwAv \y-z\ 



--.£ J^i + S J32 



with £ Jn^im'^} being first term and £ J 02(^1'^) is the remaining term of the R.H.S. For J^-^ 
using divergence theorem on dcj^izs) gives, 

^'Jl,im=[ ['[ f'^l^^±^dymliys)dy^daiz,)dysdzs 

J Me Jo JdcUeJo \y - Z\ 

dZr 



„ 'V^ ^iin%3)-n(zp). 

Oi. ''II 

J(t)£ J( 

'1 ^1 ^ ' m'(z3) 



= [ ( [ dymliys)dy^dysdzs[ V^ • ( ^ ) «., 

JiOeJo Jo Jb)c ^ \y-Z\ ' 

f fW ^.y— ^. ^ m'(z3)-iyp-Zp) 
= ^3"^3(3'3) ^dydz. 

JmJo Jo U {|3,^_2p|2+(^3_^3)2|3/2 



(A. 10) 



Setting q = {z3- ys) and dz3 = dq gives, 

1^2 



J|,(mO|<sup|m%3)|[ [ fj^^m 1(^3)1 ['"' ' '"''"/" '3/2 ^g- 

23 JcJi;J(DeJO J-ys {lyp-Zpr + g^} 



Using equation (|A7) on the inner integral gives 



Jl^{m')\<2sup\m'{z3)\ [ [ [ ^^^^ 

23 Joj£ Jai^Jo \yp Zp\ 

= 2suv\m''{z3)\\ [\dlmliy3)\dy3\\ [ [ ^f^ 

23 (Jo J I J&j£ Jwe l^p -2pl 



Z)6£^ sup|m%3)| ||a^in"||^2(o,i)^^7£^(||m"g2(n)+||53m"||j2(n)), 



„3 , 

Z3 



using equations (|A8), (|A9) and {A..2). Also we estimate JnAm^) as 



e2|j|2N (1^3(0)1 + 1^3(1)1) [ [ [ 

Joif Jo Joj 



32^ 
l^3'^3(3'3)l I53m3(3'3) 



yp-Zp\^ + yl ^J\yp-Zp\'^ + {l-y3)'^ 
<4sup|m%3)|{ C\dlml{y3)\dy3\\ [ [ p^i^l 

Z3 I Jo J I JmeJme \yp~Zp\ J 

< D8£^ sup |m"(z3)| II dim' ||^2(o,i) < ^9^^ ( || m' j^^^^ + \\ dsm' 1^2(0)) 



Z3 
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again using equations l |A.8| l, ( |A.9t and ( |A.2t . Combining estimates for J^^im^) and c/ggC/n*^) 
we get our result. D 

Recalling JgC"* '^) from eqn. ( |A.6> we write J2 = '^21 "*" '^22 "*" '^23 "*" '^24 where, 



£ t/21 



£ '^22 



Jdii)e Jo JdMs Jo 

= / / 

J 0)^(0) J a 



{\yp-Zp\^ + {y3-Z3)^}^'^' 
m'{0)-n{yp)m'iO)-n(zp) f f in' (1) ■ iiiy ^) in' (1) ■ iiiz p) 



/ / 



^ ^23-2 



W(0)-'we(0) \yp-Zp\ JwE(i)J«f(i) 

in"(0)-n(y„)in"(l)-n(zp) 



l^p-^pl 



p2 7-f 



r r^r rm''(3/3)-ra(3')'n^(0)-re(Zp) 

Jdii)F Jo [ Jw, 



and, 



^«'>{l3'p-^pP+3'|} 



,2 1 1/2 



/ 

Jb), 



ni'{y3)-n(y)m'{l)-n(zp) 

■(I){|yp-Zp|2+(1-J3)2} 



1/2 



Noting that \ni'{t)-n{z)\ = \m'^{t)\ < sup^g \ni'iz3)\ for i = and i = 1 we have using equa- 
tions (JO) and (|AJ9}, 

eV|2(mO-sup|m%3)n , " _ "^ <Dio£'(||in^||'2(n)+ ll^sm^l^a, J. (A.11) 

23 J(i>cJ(i)s\y p ^p\ 



e^Jl^iin') 



='! f! 

Jdoip JO Joif 



-m ^(0) m "(3/3) • nCy) m ^(D m "(3/3) ■ niy) 



{\yp-zp\^ + yl]"^ {\y^-zp\^ + {l-y3f]"^ 



. (A. 12) 



Note ( l^p - Zp |2 + 1 ) ^ < 1. Then eqn. ( |A.9> gives 



e^^CmO 



lla £||2 1 



<Z)ll£^(||m^||j2(a) + 



(A. 13) 



Proposition A.3. 



\j^m')\<Di2£[\\Tn'\l2(^^ + \\d3m'\l2(n))- 



Proof. As for term the J^-^ in Proposition A.2 first using divergence theorem in J^^^ from 
( |A.12t on 5a>£(3/3) we get 



eV2^4(mO 



Jo}f JoIf jo 



If m%3)-(3'p-2p)m^(0) in'iyi)-iyp-Zp)ml(l) 



{\yp-Zp\^ + yl} 



21 3/2 



{lyp-^pi'+a-ys)^} 



3/2 
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Then using m^iy^) ■ (y^ - Zp) < \yp - Zp\ supy^ Im'^Cys)] and (jA^Tf we get, 

■1 \yp-Zp\dy3 /■! \yp-Zp\d(l-y3) 



JS JoJrJoJrUOU., _« |2 , .,2)2 Jo, 



'' •'"^•'"^'•'"{lyp-ZpP + yi}^ •'°{\yp-zp\^ + {l-ys)^}'^ 

ya JiosJoii: \yp Zp\ 

and eqns. ( |A.9^ and ( |A.8^ in the last step. D 



We will now show that J2i{ni^) is the largest term in the magnetostatic terms. It contributes 
energy of 0(1) which appears in the first limit problem J^o- We split Jgi^'^'^) ^^ follows: 

'1 m''{y3)-n{y ) m%Z3)-n{Zp) 



I fl f 

Jdo)f: JO Jdoip JO 



^^V'"')-;... ;„ L L ^y.^. 



r r^ r r^- m'{y3}-n(yp)m''{y3}-n(zp) 

JdoJeJO Jdoj.Jo £^\y-z\ 

■1 f f'^m%y3)-n{yp){m''iy3)-m'iz3))-n{Zp) 



■/ II I 

JdoipJO JduicJO 



£'^\y-z\ 



= Jn-,Am'} + JLoim'). 



Next we show the following proposition. 
Proposition A.4. 

k212(*"')| ^Di3£^'^ [ ||m^||j2(n)+ 11^3^^11^2(0)). 

Proof. Using Divergence theorem in y, variable as in ( |A.10^ and Fubini's theorem we get, 

/. /.I /.I /■ m'^{y3)-n(y„) 
^%'i2= / {m'{z3)-m%y3))-n{zp)da{yp) 

Jaa,,J0 Jo Ja^,^|3,^_;jp|2 + (3,3. 23)2 

r fl fl r m'{y3)-(zp-yp) 
= / / / / ]- ^- — ^(m%3)-m (j3))-n(zp)dy^ 

Jda>,JO Jo JoJe {\yp-Zp\^ + iy3-Z3)^ \ 

f f ~r ^ f'-^r ^r , f^ in^(z3)-m%3) 

Jaw,J(«, Jo Jo ||y^_2;^|2 + (^3_23)2 I 

Now note that — -^g — ^ ^'^'^ |w(Zp)| = 1 which gives 

{lyp-2pl^+(3'3-Z3)^ } 

.2 7£ .^£... ^^l l^£... xl If f f^ f^ \m%Z3) -m%y3) 

■iy3-z3)^ } 

(A. 14) 



'Jl,,im%y3))\<snp\m^iy3)\ \[ [ [ [ l^^("3) ■- ----■^^ ^,3^. 
js Usw^JweJo Jo {l3'p-Zpr + 
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Note that, 

1 11 
< (j\ ]^5) 

{\yp-Zp\^ + {y^-Zi)''} l3'p-Zp|l/4|3;3-23|V/4 

Then 

ri r Mm^(23)-m%3)lrfz3rfy3 l r^ r^ \m'(zs)-m'{yi)\ 

Jo Jo {\yp-zp\' + iy3-z3)'} ~\yp-Zp\^^Jo Jo lys-zal''^ '' ^' 

^ 1^ _\ |l/4 ll^3^'(y3)Ll(0,l) (A.16) 

because of the fact that I I "^ ^^_ 7^/4^^ dz^dys denotes the seminorm in the fractional 

Sobolev space W4'^(0, 1) and by the continuous embedding of W^'^iO, 1) c Wi'^iO, 1). The 
integral above cannot be bounded by norm in W ' alone unless m^ is a constant, which 
is shown by the surprising result Proposition 1 in jBrezis, 2002 J . Also note using Holder 
inequality 

II^3'"1Il1(0,1) = / \d3m'\X(0,l)dy3 < ||I(0,1)L2(o,i)||^3'"1Il2(0,1) = II^3'"1Il2(0,1)- 

Then using ( |A.16t in eqn. ( |A.14| I along with the above result we get. 



je r--er ^^^ l-^r mU f f f \m'(z3)-m'(y3)\ ^ 

J2i2(m (3/3)) < sup m (3/3) < Y. ;^- Wjdzs 

ys [Jda,,Jco,Jo Jo {lyp-Zpr + iys-zaV \ 

l-^r M ll;,y-^ll i f f dypdaizp} ^ 

<SUp|m (J3)| . ia> 1^2(0,1) 1^ _^ ,1/4 

ys ijawfjwc \yp Zp\ J 

= Z)i4 £^ e^/^ ( \\m'\\l2^^^+\\dsm'\\l2^a)) 

using calculation like in eqn. jA.S) to get the £ £ term and eqns. ( |A.9| I and (|AJ2). D 



In 2-dimensional micromagnetics on a domain Y £ IR for a constant magnetization m £ 
i? (T, TOsS ), the demagnetization field is given by, 

r X ~ y 

hm(x)= f -^m-n{y)dy (A.17) 

Jdw \x-yr 

and magnetostatic energy is given by, 

S2d= [ [ m- ^~^^^ mn{y)dy. (A.18) 
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Proposition A.5. 

J2ii{m')-2ji\a)i:\ I \ml{y3)fdy3 <Z)i5 e(||in''||^2(n)+ ||^3'ra''||i2(^))- 
Proof. Using the Divergence theorem on Zp as in ( |A.10| I and a subsequent change of variables 



qiza) = Z3- ya, followed by ( |A.7| f (as in Proposition |A^ we get 

'1 r^ rm'{y3)-n{yp)m^{y3)-n(zp) f f^ /"i fm'{y3)-n{yp)m^{ys)-{yp-Zp) 



I Hi- 

Jdoip JO JO Jdoif- 



^p - Zp P + (3/3 - Z3)2 



/ III 

JdiOp JO JO Jch 



{|yp-Zp|2+ (3/3-23)2} 



3/2 



= fff m'iy3)-n{yp)[ 

Jdhif: JO JcOe J- 



i-j-s m%y3)-{y„-Zp) 



^3 {\yp-Zp\^ + q^} 



3/2 



dq 



1 r m''{y3)-iyp-Zp) t y3m''iy3)-n{yp) (1- y3) m''iy3)-n{yp) 



r r^ r m (^3) 

Jdii)i: Jo Jii>£ I yp ■ 



yl + \yp-Zp\^ Jil-ysY + lyp-Zpl 



= :£2[J2'lll(^') + ^2112(^')) 

where e^J^^-^-^im") = f If 

Jdo)^ Jo Jih 

ing term. Set J^{fn^) as 



1 r mHysHyp-Zp) ys mHy3)n(yp) , 2 re /-■— -£\ 4.1. 

' — ; ^ — , =^^ and £ J L -.Am ) the remain- 



jm') 



f^ f f —e, ^ ~, m'{y3}-{yp-zp} 
'■=11 I "I iy3)-n{y ) ■ -^ daiy )dzpdy3, 



\yp-zp\ 



< < 



(A. 19) 



Let i? := max 2£ ^\Xp-yp\, [zp e a)E,yp e dcUf-), andnote 

' \yp-Zp\^ r- ^ D 

^^, {ory3>R£, 
1 for 3/3 <i?£ . 

Noting that |re| = 1, \m^{y3)-n{yp)\ < |m''(3'3)| and \m^{y3)-(yp-Zp)\ < \m'{y3)\ \yp-Zp\ 

ys 



ya 

yl + \yp-zp\^ 



'^m'{y3)-(yp-Zp)^^ 



rm (3/3 j 

JO \yp 



-m 



p I 
mHy3) 



(y3)-n{yp)[^ 



\lyl+\y. 



dy3 



p ■^p' 



p ^p 



yl + \yp-zp\^ 



< r 1^ 

~ Jo \ y 

^ I ] :dy3 + / m (3/3)! 

Jo \ yp-Zp\ Jrc \ 



^ \in'{y3)f 



r m 

jRe I y 



ly^-gpi 
2y? 



p ^p 



dy3 



ys 



yi + \yp-zp\^ 



< sup |m^(3'3)|' 



I Jo \yp-Zp\ jRe lyl ] 
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Using above result and noting that ^3 [y^ ) = -y^ we get 



8^ 



^mHy3)-iyp-Zp) _ r 3/3 



\ r r r^m \y3)-\y„-Zr,) , 

•'«-•'-'' Alii ^;rf^«'<»>-«^'>( 



' -'• ■ i/yf + iy,-',!"' 



yz J dole J<^e ( Jo I ^p " ^p I JRe ^^3 

_,^ ^ n r r 1 R£ \yp-zp\ 1 1^ 

<sup|m%3)r \- T-^ 

3-3 JdciteJoJe U^p Sip | -^ 



:y3 JdojEJoE ^\yp-Zp\ Z ZK£ I 

Note from equation jA.S) , the term | dzp = -Dise. So the first integral above is 

Jtt)^ \yp ~ Zp\ 

Re I I ■ -dzp ~ DiqE^. The second integrand is 0(e) and so its integral is of 0{£^). 

Jd(U£ J(i>E \yp - Zp\ 

The third integrand is bounded by 1, since by definition i?e > ly^ - Zp |. So the third integral 
I f ^ ^ dzp ~Dn£^.So\Jl-J2iii\^Dif,£svi^y^\m'^(y3)\'^.J2ii2C&n'beireaieAi\ie 



J2ii('n'')-2Jo(in'') |<-Di8e sup|m''(3/3)l^ <-Di9£(||'n'^||i2(Q) + ||33m''||^2(n))- 



same way to give the result on using eqn ( |A.2t 

>\m-(y3)\- su-ig£(\\m-{ ' "" "" 

J'3 

We get our result noting that JQ{m) is exactly the 2-D magnetostatic energy S2d defined in 
( |A.18^ and for a circular cross-section w^ it is well know that 



S2dim'') = Jo{m'') = 7i\a)e\j liripiya)] dys = £^Ji\a)\ j \mp{x3)\ dx^. 



n 



A.3 Final Estimate for S^im') 

The excange energy of ?n'^ is bounded by equation ^5.3) , -j- > £~^||Vp/n^||^2m) + ||53in'^||^2(Q)- 
Using Remark [A.l| we get first, 

Combining Propositions |A.l| , |A.2| , |A.3| , |A.4| , |A.5| and and equations l |A.llt and ( |A.13| I we get, 

<g'^{m')-7t\a)\f \mliys)fdys = {Oi£) + Oi£^'^))[\\m'\\l2^n^ + \\d3m'\\l2^a^]. 
Combining the two we get 

(f^dnO -n\a,\f Imliyafdya = 0(e) + 0(£^'^). (A.20) 
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Remark A.2. TTie above result ( |A.20| I is true for any magnetization m as long as the magne- 

-9 11 l|2 II ||2 

tization satisfies the exchange bound K^ > £ \\'^pfn\\i,2im + II^S'^llL^fn)- 

Let Tn° be a constant vector in nisS . If /re" is the rescaled version of m", recall the result 
in equation ( |3.6t gives, 

S^,im'') = e-'E,,^,, = n'\mlf-e^(\mlf-2\mlf] + 7:h'(^-\mlf] 



where we define Qi and Q2 as in equation ( |3.6> . 

Proposition A.6. Let m° be a constant on rngS'^ and H^{Q.E;msS^)3Tn'^ = ni° +M . Then 
the following holds in terms of the rescaled magnetizations [Tn^,Tn",M'^], 

|2 



t J LI 

Proof First note ^gjO as given in ( |3.6t on rescaling to ^^o gives 



?l7r 


ire° 


16e 

H 


■ 'n; 


+ nf^ 


K ■ 







3 


[-'Zml\ 




[-^m^g] 



Note SI ptn^ = Vp(ire° +M'^) = SI pM'^ as ire" is constant. Lemma A.l gives along with Young's 
inequality gives 



^^^PpM'tnny 



^0^' 11,^. II 2 



Using propositions |A.l| , |A.2| , |A.3| , |A.4| , |A.5| and equations ( |A.llt and ( |A.13t we get 



-D,,£^'\\M^t,^^^ + \\d,M^t,^^^. 



Adding the two togethc 


jr we get, 




d 


n 


+ <§'^(M")- 


-^/o^^^ 








^' L2(n,-^ 



1)20£^'^(|| 



^1li2(n)+||53M^|li2(„))- (A.21) 
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Note by the linearity of Maxwell's equation 

Dividing by £^, noting that S^im^) = ^ ^ /^s {hj^s \ dy, and rescaling m^ & M , 

- £^7i I (m° • Mp - 2mlM'^)dx (A.22) 

where we have used Young's inequality to bound the last two terms in the final step i.e. 

D2i£\ml-Ml\ < {D2i£^\ml\). (VA |M;|) < ^\m°/ + ^\Mlf and similarly for 
the TOqAfo term. Also n I |m!,| = n I |m°| +n( |M!,| +2nl m" • Ml, . Subtracting the 
two and noting ||'w°||^2(q) = mg\D.\ since Im"! = mg, we get 

g'^(m')-n[ \mlf>£Qi + £^Q2 + S'd(M')-7i[ \M'/ - j\\M'\\l2^^^-D2s£^m^,\n\. 

J LI ^ LI ^ 

Using this and eqn. ( |A.21> we get ( note Vpin" = VpM"^ ) 

c J LI 

>.Qi + .^Q2 + ^||VpM^||j2,„, + ^(||M^||j2(„, + ||a3M^||j2(„,)-Z)i8£2 
for £ small enough. D 



Remark A.3. If m" = (0,0, nis), we get a simpler estimate than in above Proposition A. 6 



Note m' = m''+M' gives \m'f = m^ = \m°\^ + \M' \^ + 2m'' ■ M' = mf + \M'f + 2mlMl 
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which means -2m'^M^ = \M'^\^ and nip -Mp = 0. Substituting these in ( |A.22t we get 

-£^n [ {m''-Ml-2mlMl)dx 
Jn 

= S',{m'') + S^{M') + 2ji [ m'i-Mldx+[—-£^n) [ iM'l^dx. 

Jn 3 Jn 



Using the above and proceeding with the remaining part of the estimate in Proposition \A. 6\ we 
get the following result, 



J-'\' 



t- J LI 

> eQ, + e^Q2 + ^ PpM'\\l2,n, + f ( Il^1li2(^) + \\d3M' Wlnm) (A.23) 

for £ small enough. 
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